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ABSTRACT 

The  dependence  of  the  hole  density  on  that  of  Mn  sites  in  Gai_xMnxAs  is  studied  within  a 
mean-field  approach  to  the  hole-mediated  ferromagnetism  in  III-V  Mn-based  semiconductor 
compounds.  We  parametrize  the  hole  concentration,  p,  as  a function  of  the  fraction  of  Mn  sites,  x, 
in  terms  of  the  product  m*  (Jpd)2  (where  m*  is  the  hole  effective  mass  and  J|X|  is  the  Kondo-like 
hole/local-moment  coupling),  and  the  critical  temperature  Tc.  By  fitting  m (Jpd)'  to  experimental 
data  for  Tc(x),  we  establish  the  dependence  of  p on  x,  which  is  interpreted  in  terms  of  a reentrant 
metal-insulator  transition  taking  place  in  the  hole  gas. 


INTRODUCTION 

The  discovery  in  the  early  1990's  of  ferromagnetism  in  III-V  materials  alloyed  with  transition 
elements  like  Mn  [1,2]  has  increased  the  interest  in  the  study  of  the  electronic,  optical  and 
transport  properties  of  diluted  magnetic  semiconductors  (DMS).  Ferromagnetic  semiconductors 
bring  about  the  possibility  of  controlling  both  spin  and  charge  degrees  of  freedom.  The 
combination  of  this  feature  with  the  capability  of  growing  low-dimensional  structures  has  opened 
up  new  perspectives  for  the  production  of  spintronic  devices.  Applications  would  include  non- 
volatile memory  systems  [3]  and  quantum  computing  [4]. 

A great  deal  of  attention  has  been  focused  on  Gai.xMnxAs  alloys,  which  exhibit  very 
interesting  magnetic  and  transport  properties.  When  Mn  atoms  replace  Ga  in  GaAs,  their  five  3d 
electrons  remain  localized  in  a core  state,  giving  rise  to  S = 5/2  local  moments.  In  addition,  Mn 
atoms  have  one  less  electron  in  the  4s  level  than  Ga,  so  they  act  as  acceptors  generating  hole 
states  in  the  material;  for  large  enough  doping,  these  states  merge  to  form  an  impurity  band.  It 
should  be  pointed  out  that  the  equilibrium  solubility  of  Mn  atoms  in  GaAs  is  quite  low  [5],  being 
only  of  the  order  of  10!v  cm'3.  However,  with  help  of  molecular-beam  epitaxy  techniques  at  low 
temperatures,  it  has  been  possible  to  produce  homogeneous  samples  of  Gai_xMnxAs  with  x as  high 
as  0.071.  Magnetization  measurements  in  these  systems  show  that  for  0.015  < x < 0.071  they 
become  ferromagnetic,  with  doping-dependent  critical  temperatures  Tc(x)  reaching  a maximum  of 
1 10  K for  x = 0.053  [6], 

The  origin  of  ferromagnetic  order  in  these  materials  is  not  fully  understood,  but  there  exists 
consensus  on  the  fact  that  it  results  from  the  exchange  coupling  between  the  localized  Mn 
moments  mediated  by  the  holes.  The  strength  of  the  of  the  coupling  is  expected  to  depend  on  the 
hole  concentration  p,  which  eventually  determines  the  Curie  temperature  Tc(x).  In  principle,  one 
would  expect  that  each  Mn  would  provide  one  hole,  leading  to  a density  of  holes  equal  to  that  of 
the  magnetic  ions.  However,  this  has  not  been  confirmed  by  experiment.  Although  an  accurate 
determination  of  the  hole  concentration  is  hindered  by  the  anomalous  Hall  term,  experimental  data 
indicate  that  p is  of  the  order  of  a 15  to  30  % fraction  of  that  of  magnetic  ions  [6-9].  The 
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mechanism  responsible  for  the  discrepancy  between  hole  and  Mn  densities  is  not  clear.  Matsukura 
et  al  [7]  have  pointed  out  that  such  discrepancy  might  be  due  to  compensation  of  Mn  acceptors  by 
deep  donors  such  as  As  antisites,  which  are  known  to  be  present  at  high  concentration  in  low- 
temperature  grown  GaAs  [10].  However,  other  mechanisms  such  as  the  formation  of  sixfold- 
coordinated  centers  with  As  (MnfiAs),  which  would  compensate  Mn  atoms  on  substitutional  Ga 
lattice  sites  [11],  should  not  be  ruled  out.  As  a consequence,  no  theory  relating  the  hole 
concentration  and  that  of  Mn  has  as  yet  been  proposed.  However,  this  issue  is  of  great  interest  for 
the  design  of  new  devices.  The  main  purpose  of  the  present  work  is  to  present  a preliminary 
quantitative  analysis  on  this  matter. 


THEORETICAL  FRAMEWORK 


The  appearance  of  ferromagnetism  in  Gai.xMnxAs  has  been  explained  on  the  basis  of  the 
following  picture.  The  Mn  ions  interact  with  the  holes  via  a local  antiferromagnetic  Kondo-like 
exchange  coupling  Jpj  between  their  moments  [9,12-15],  which  results  in  the  polarization  of  the 
hole  subsystem;  this,  in  turn,  gives  rise  to  an  effective  ferromagnetic  coupling  between  the  Mn 
moments.  We  point  out  that  there  has  been  some  debate  as  far  as  the  details  of  the  above  picture 
are  concerned  (e.g.,  whether  or  not  such  effective  interaction  is  well  described  by  an  RKKY  term 
[16,17]).  Nevertheless,  there  is  consensus  on  the  fundamental  role  played  by  the  hole-mediated 
mechanism.  We  emphasize,  however,  that  the  approach  we  follow  here  does  not  depend  on  the 
details  of  the  effective  Mn-Mn  interaction. 

We  can  thus  think  of  two  coupled  subsystems  (localized  moments  and  holes)  as  described  by 
the  total  Hamiltonian 


H = HMn  + H„+JpdXS,.s,5(ri-RI).  (1) 

u 

Here  HM„  describes  the  direct  (i.e.,  non-hole-mediated)  antiferromagnetic  exchange  between  Mn 
spins,  Hh  describes  the  hole  subsystem,  and  the  last  term  corresponds  to  the  aforementioned  Mn- 
hole  exchange  interaction,  with  Si  and  Si  labeling  the  localized  Mn  spins  (S  = 5/2)  and  the  hole 
spins  (s  = 1/2),  respectively.  As  a first  approach,  we  neglect  HMn  and  consider  Hh  within  a 
parabolic-band  effective-mass  approximation;  a more  general  description  of  Htl  may  include  a 
multi-band  Kohn-Luttingcr  treatment  [9, 1 8],  effects  of  impurity  potentials,  a site  energy  term 
arising  from  the  Mn  potential,  and  a correlation  energy  representing  hole-hole  repulsion. 

The  problem  defined  by  the  Hamiltonian  in  Eq.  (I)  is  indeed  quite  complex,  but  some 
physical  insight  can  be  gained  by  a mean-field  approach.  Assuming  that  the  magnetization  density 
of  the  hole  subsystem,  Mh  = < nr  - ni  >,  is  uniform  within  the  length  scale  of  magnetic 
interactions,  the  Mn  magnetization  is  given  by 


M = NM„gpBM,  = nMllgpnS  Bs 


f JP,S  ' 
2k  bT 


M, 


(2) 


where  nMll  = x n,  is  the  density  of  Mn  ions,  with  n,  being  the  density  of  Ga  lattice  sites,  Mi  is  the 
magnetization  density  of  the  Mn  ions,  g = 2 is  the  Mn  Lande  g - factor,  and  Bs  is  the  Brillouin 
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function.  The  assumption  of  homogeneity  of  the  hole  magnetization  allows  us  to  evaluate  Mh  self- 
consistently  by  considering  a Fermi  sea  of  holes  with  effective  mass  m‘,  in  the  presence  of  the 
mean  magnetic  field  generated  by  the  Mn  ions.  We  therefore  find  that 

Mh  °=  m‘'Jpd  x M,  pl/3.  (3) 


RESULTS  AND  DISCUSSION 

The  critical  temperature  as  a function  of  hole  density  and  Mn  composition  is  obtained  by 
linearizing  the  self-consistency  relations  given  by  Eqs.  (2)  and  (3),  i.e., 


which  leads  one  to  write  the  hole  concentration  as 


P = ?- 


Tc(x) 


(m  */mc)Jpdx 


(5) 


where  £ = 5.29  x 10  ls,  in  units  such  that  Jpd  is  given  in  eV  nm3. 

The  sought-after  relation  between  p and  x is  provided  by  Eq.  (5).  However,  it  depends  on  the 
parameters  m*  and  Jpd,  whose  values  have  not  been  determined  with  certainty  from  the 
experimental  data.  If  one  tentatively  uses  m*  = me,  J = 0.15  eVnm3  [14],  and  the  experimental 
transport  data  [6]  for  Tc , the  resulting  estimates  for  p(x)  are  as  given  by  the  full  circles  in  Fig.  1 . 
The  error  bars  reflect  the  uncertainties  in  the  determination  of  Tc,  as  displayed  in  Fig.  3(c)  of  Ref. 
6.  Nevertheless,  a more  reliable  estimate  can  be  obtained  by  adopting  the  following  strategy. 

Since  the  Hall  resistance  measurements  [19]  yield  an  unambiguous  [20]  value  of  p = 3.5  x 102n 
cm"3  for  the  sample  with  x = 0.053,  for  which  Tc  = 1 10  K,  one  may  use  Eq.  (5)  to  determine  the 
value  of  the  product  (m‘/me)  (JPd)2  = 2.4  x 10'2  (eV  nm3)2.  Thus,  by  taking  this  estimated  value 
along  with  the  experimental  data  [6]  for  Tc  into  Eq.  (5),  one  gets  the  data  for  p(x)  shown  as  full 
squares  over  a wide  range  of  x.  The  adequacy  of  the  above  procedure  is  confirmed  by  the  fact  that 
the  calculated  values  for  p(x)  lie  below  the  concentration  of  Mn  ions,  shown  as  a dotted  line,  in 
agreement  with  experiment.  Before  accepting  these  estimates  for  p(x)  at  face  value,  one  should 
note  that  a closer  look  at  the  experimental  data  for  Tc(x)  [6]  suggests  a linear  behavior  in  the 
range  of  x of  the  order  0.015  - 0.035,  which  would  imply,  through  Eq.  (4),  a constant  p in  that 
range;  such  constant  behavior,  however,  should  not  prevail  at  low  concentrations,  x — > 0,  where 
presumably  one  should  have  p <*=  x — > 0.  These  considerations  are  incorporated  in  the  full  curves 
displayed  in  Fig.  1 , which  lie  within  the  corresponding  error  bars  of  the  calculated  hole 
concentrations. 

The  boundaries  of  the  metal-insulator  transitions  (MIT’s),  as  determined  from  resistivity 
measurements  [6],  are  indicated  in  Fig.  1 by  vertical  dashed  lines.  We  remark  that  the  present 
theoretical  estimates  for  p in  the  insulating  phases  are  based  on  the  assumption  that  the  scale 
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Figure  1.  Theoretical  results  for  the  hole  concentration  as  a function  of  the  fraction  of  Mn  sites 
for  (Ga,Mn)As  ferromagnetic  alloys.  The  dashed  line  corresponds  to  a hole  concentration  equal  to 
that  of  Mn  sites,  whereas  the  filled  circles  and  squares  are  the  present  mean-field  results.  Full 
curves  are  obtained  as  reasoned  in  the  text.  I and  M respectively  denote  insulating  and  metallic 
phases. 

of  localization  length  in  insulating  samples,  though  finite,  is  significantly  larger  than  the  length  of 
magnetic  interactions  [9],  in  which  case  the  present  mean-field  approach  is  a good  starting  point. 

The  most  prominent  feature  of  Fig.  1 is  the  fact  that  p(x)  displays  a maximum  within  the 
metallic  phase.  It  means  that  attempts  to  increase  Tc  should  be  carried  out  for  samples  in  the 
metallic  phase,  for  Mn  concentrations  about  0.05.  Moreover,  notwithstanding  the  considerable 
uncertainties  [21]  in  the  measurements  of  p,  the  present  theoretical  estimates  for  the  hole 
concentration  exhibit  trends  in  quite  good  agreement  with  those  obtained  from  Hall  measurements 
by  Matsukura  et  al  [7]. 

The  fact  that  the  dependence  of  p with  x is  essentially  related  to  the  occurrence  of  MIT's 
taking  place  in  the  hole  subsystem  can  be  seen  from  the  following  qualitative  picture.  Within  our 
approximation,  the  Fermi  energy  tracks  the  behavior  of  p,  since  £r  a p2/\  while  the  exchange 
splitting  A a x.  Figure  2(a)  shows  the  schematic  impurity  bands  for  each  spin  channel,  in  the  very- 
low  doping  regime  in  which  the  gas  is  supposed  to  be  unpolarized.  As  x increases,  the  gas  can 
sustain  polarization  and  still  be  insulating,  provided  the  Fermi  energy  lies  below  the  mobility  edge, 
as  shown  in  Fig.  2(b).  Further  increase  in  x causes  eP  to  increase  and  to  lie  within  the  delocalized 
states  of  the  up-spin  impurity  band,  as  depicted  in  Fig.  2(c):  The  system  becomes  metallic. 
Whether  or  not  the  Fermi  energy  also  lies  within  the  delocalized  region  of  the  down-spin  impurity 
band  is  a very  interesting  question,  which  cannot  be  answered  by  our  simple  model;  the  solution 
of  this  particular  issue  should  have  bearings  on  the  efficiency  of  Ga,  .x  Mn,As  - based  devices  as 
spin  filters.  Once  en  reaches  a maximum  within  the  metallic  phase,  its  initial  decrease  upon 
increasing  x is  compensated  by  an  increase  in  A.  so  that  the  Fermi  level  still  lies  within  the 
delocalized  states.  However,  with  continuing  increase  in  x the  exchange  splitting  can  no  longer 
make  up  for  the  decrease  in  £K,  and  the  latter  eventually  crosses  the  mobility  edge  again,  lying 
within  localized  states  [Fig.  2(d)]:  The  system  reenters  an  insulating  phase. 

The  simple  approach  used  in  the  present  work  provides  a straightforward  and  immediate 
estimate  of  the  hole  concentration  in  Ga,.,  Mn,As  as  a function  of  both  x and  Tc  [22],  which  is  an 
important  ingredient  to  guide  further  experimental  work  on  this  subject.  Clearly,  the  present 
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approach  could  be  extended  to  include  a more  complete  description  of  the  acceptor  states,  taking 
into  account  the  spin  degrees  of  freedom,  spin-orbit  coupling,  compressive/tensile  strains,  etc. 
Moreover,  a proper  treatment  of  disorder  — e.g.,  by  explicitly  considering  a random,  instead  of 
continuous,  distribution  of  Mn  ions  — should  lead  to  a more  realistic  description  of  the  MIT.  In 
this  respect,  many-body  effects  due  to  correlation  among  the  holes  should  also  influence  p(x), 
especially  in  the  insulating  phase.  Of  course,  an  appropriate  description  of  the  physical 
mechanisms  related  to  As  antisites  and  Mn';As  centers,  as  far  as  the  hole  vs.  Mn  concentrations  is 
concerned,  is  certainly  a formidable  task,  which  nonetheless  deserves  future  theoretical  attention. 

In  summary,  we  have  considered  a mean-field  theory  for  the  hole-mediated  coupling  between 
the  localized  Mn  moments  in  Gai.xMnxAs.  The  hole  concentration  as  a function  of  the  fraction  of 
Mn  sites  is  obtained  in  terms  of  the  product  m*(Jpd)2  and  the  critical  temperature  Tc.  By  using 
experimental  data  for  these  latter  quantities,  we  have  established  that  the  occurrence  of  a reentrant 
MIT  taking  place  in  the  hole  gas  can  be  traced  back  to  the  dependence  of  the  hole  concentration 
with  x.  Although  not  shown  here,  we  have  also  calculated  the  dependence  of  the  Mn 
magnetization  with  x,  for  different  temperatures,  and  found  that  as  T increases,  the  width  of  the 
composition-dependent  magnetization  decreases  drammatically,  and  that  the  magnetization 
maxima  also  decreases,  indicating  the  need  for  quality-control  of  Mn  doping  composition  in  DMS 
devices. 
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Figure  2.  Schematic  density  of  states  (DOS)  versus  energy  for  the  impurity  band,  for  up-  (top) 
and  down-spins  (bottom).  Under  each  DOS  curve,  the  hashed  and  empty  regions  correspond, 
respectively,  to  delocalized  and  localized  states;  these  are  separated  by  mobility  edges.  The 
exchange  splitting  is  proportional  to  the  off-set  between  the  T and  4-  bands,  and  the  Fermi  energy 
(ef)  increases  to  the  right. 
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